Fourth-order algorithms for solving the multivariable Langevin equation and the Kramers equation.

نویسندگان

  • H A Forbert
  • S A Chin
چکیده

We develop a fourth-order simulation algorithm for solving the stochastic Langevin equation. The method consists of identifying solvable operators in the Fokker-Planck equation, factorizing the evolution operator for small time steps to fourth order, and implementing the factorization process numerically. A key contribution of this paper is to show how certain double commutators in the factorization process can be simulated in practice. The method is general, applicable to the multivariable case, and systematic, with known procedures for doing fourth-order factorizations. The fourth-order convergence of the resulting algorithm allowed very large time steps to be used. In simulating the Brownian dynamics of 121 Yukawa particles in two dimensions, the converged result of a first-order algorithm can be obtained by using time steps 50 times as large. To further demonstrate the versatility of our method, we derive two new classes of fourth-order algorithms for solving the simpler Kramers equation without requiring the derivative of the force. The convergence of many fourth-order algorithms for solving this equation are compared.

برای دانلود متن کامل این مقاله و بیش از 32 میلیون مقاله دیگر ابتدا ثبت نام کنید

ثبت نام

اگر عضو سایت هستید لطفا وارد حساب کاربری خود شوید

منابع مشابه

0 Ju n 20 00 Fourth Order Algorithms for Solving the Multivariable

We develop a fourth order simulation algorithm for solving the stochastic Langevin equation. The method consists of identifying solvable operators in the Fokker-Planck equation, factorizing the evolution operator for small time steps to fourth order and implementing the factorization process numerically. A key contribution of this work is to show how certain double commutators in the factorizat...

متن کامل

The smoothed particle hydrodynamics method for solving generalized variable coefficient Schrodinger equation and Schrodinger-Boussinesq system

A meshless numerical technique is proposed for solving the generalized variable coefficient Schrodinger equation and Schrodinger-Boussinesq system with electromagnetic fields. The employed meshless technique is based on a generalized smoothed particle hydrodynamics (SPH) approach. The spatial direction has been discretized with the generalized SPH technique. Thus, we obtain a system of ordinary...

متن کامل

The Langevin or Kramers Approach to Biological Modeling

In the Langevin or Ornstein-Uhlenbeck approach to diffusion, stochastic increments are applied to the velocity rather than to the space variable. The density of this process satisfies a linear partial differential equation of the general form of a transport equation which is hyperbolic with respect to the space variable but parabolic with respect to the velocity variable, the Klein-Kramers or s...

متن کامل

A Note on the Smoluchowski-kramers Approximation for the Langevin Equation with Reflection∗

According to the Smoluchowski-Kramers approximation, the solution of the equation μq̈ t = b(q μ t ) − q̇ t + Σ(q t )Ẇt, q 0 = q, q̇ 0 = p converges to the solution of the equation q̇t = b(qt) + Σ(qt)Ẇt, q0 = q as μ → 0. We consider here a similar result for the Langevin process with elastic reflection on the boundary.

متن کامل

High-accuracy alternating segment explicit-implicit method for the fourth-order heat equation

Based on a group of new Saul’yev type asymmetric difference schemes constructed by author, a high-order, unconditionally stable and parallel alternating segment explicit-implicit method for the numerical solution of the fourth-order heat equation is derived in this paper. The truncation error is fourth-order in space, which is much more accurate than the known alternating segment explicit-impli...

متن کامل

ذخیره در منابع من


  با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید

برای دانلود متن کامل این مقاله و بیش از 32 میلیون مقاله دیگر ابتدا ثبت نام کنید

ثبت نام

اگر عضو سایت هستید لطفا وارد حساب کاربری خود شوید

عنوان ژورنال:
  • Physical review. E, Statistical, nonlinear, and soft matter physics

دوره 63 1 Pt 2  شماره 

صفحات  -

تاریخ انتشار 2001